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We make a parallel excursion to the generalised Maxwell and Kelvin-Voigt models to find which 
general conditions may lead to the weak power law behaviours of the elastic and viscous moduli, 
G'{w) ~ G"{vu) ~ w°‘, with 0 < a < 1. These occur in a wide variety of soft materials, including 
living cells. We find that these behaviours are obtained for analogous, but essentially different 
power law spectra of relaxation and retardation times, stemming from the different nature of the 
models. The first, based on the generalised Maxwell model, retains the fluid properties of a soft 
glassy material (SGM) above the glass transition. The second, based on the generalised Kelvin- 
Voigt model, retains the elastic properties of a soft solid or gelled material. Using these models, we 
calculate with accuracy the power law exponents of G'{w) ~ and G"{w) ~ , for the whole 

families of relaxation and retardation spectra previously obtained. Within the region corresponding 
to 0 < a < 1, we find x > y ~ a and y > x ~ a for respectively the generalised Maxwell and 
Kelvin-Voigt models. Beyond this region, we find very different and interesting exponents. 


INTRODUCTION 

In the last couple of years, we have studied the rheo¬ 
logical properties of Staphylococcus aureus cultures m- 
We have observed a very rich viscoelastic behaviour, both 
under steady-state and oscillatory flows, as a consequence 
of bacteria multiplication and density dependent aggre¬ 
gation. The elastic and viscous moduli, G' and G", ex¬ 
hibit approximately the same weak power-law behaviours 
with the angular frequency, w, of the oscillatory defor¬ 
mations. Interestingly, the exponents of these power-laws 
depended on the growth stage of the culture, starting to 
be close to 0.80 at the beginning of the exponential phase 
of growth, when the cell density was still relatively small, 
decreasing to 0.23 at the end of the exponential phase, 
when the cell density became large, and increasing again 
to 0.39 at the late phase of growth. 

Microrheology measurements performed by attaching 
magnetic microspheres to the cell and applying external 
oscillatory magnetic flelds also revealed the same type 
of power law behaviours [3]. The weak power law expo¬ 
nents, which are associated in this case with the mechan¬ 
ical properties of the cell cytoskeleton, are often close to 
0.2, for different types of cells [3], but may take, in some 
occasions, higher or lower values [5]. 

In fact, long before this set of experiments, it had al¬ 
ready been recognised that weak power law behaviours, 
G'{w) ^ G''{w) ^ (0 < a < 1), occur in a large va¬ 

riety of soft materials. In diluted unentangled polymers, 
the power law exponents are associated with the relax¬ 
ation modes of each molecular chain [Ml. In pre-gel and 
post-gel polymers, they are associated with percolation 
and network self-similarity nnHni ( an extensive review is 


found in [TS] ). In the context of colloidal glasses [MU, 
they are predicted by the mode coupling theory (MCT), 
and are associated with the fast relaxation modes at¬ 
tributed to the localised motion of particles inside cages 
that are formed by their neighbouring particles. The uni¬ 
versality of the weak power law behaviours, even for sys¬ 
tems devoid of any obvious fractal structure, like “foams, 
emulsions, pastes and slurries”, led Sollich and coworkers 
to propose the phenomenological “Soft Glassy Materials” 
(SGMs) model [22JI23]' Based on the idea of structural 
disorder and metastability, common to all SGMs, this 
model relates the power law exponent to a mean-field 
noise temperature x = a -\- 1, with a glass transition oc¬ 
curring aX X = 1. But the literature in this held is vast 
and constantly increasing. Other important theoretical 
models have been consecutively proposed, depending on 
the particular characteristics of the materials in analysis 
(e.g. [U)- B would be very difficult for us, and clearly 
beyond the scope of our work, to give a comprehensive 
list of all these models. 

In the context of cell mechanics, the weak power law 
behaviour of the cytoskeleton has been recently associ¬ 
ated with the “Glassy Worm Like Ghain” (GWLC) model 
[U US] , so-called from its analogy to SGM model (see 
m for a review). This model defines an average sepa¬ 
ration between the crosslinks along the filaments. If the 
relaxation modes have a wavelength shorter than this 
separation, its relaxation time follows the “Worm Like 
Chain” (WLC) model for semiflexible polymers, which 
largely compose the cytoskeleton. Otherwise, the relax¬ 
ation spectrum is stretched through an effective Boltz¬ 
mann factor with a characteristic energy e that must be 
overcome to induce a conformational change of the net¬ 
work. This model predicts a high frequency regime with 
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a power law exponent a = 3/4 (corresponding to WLC 
model) and a low frequency regime with a second power 
law exponent that depends on the phenomenological pa¬ 
rameter e. Extensions of this model deal with the possi¬ 
bility of transient crosslinking between filaments |26j . 

A large majority of the models uses a microscopic de¬ 
scription of the materials to propose a specific spectrum 
of relaxation modes. Thus, they may be directly asso¬ 
ciated with the generalised Maxwell model (see Fig. 
top) from the linear theory of viscoelasticity [551 US] • As 
it is known and will be shown also in this article, there 
is a well established connection between the weak power 
law behaviours previously mentioned and a spectrum of 
relaxation modes that behaves itself as a power law of 
its relaxation times. However, the cell cytoskeleton has 
also been previously associated with a generalised Kelvin- 
Voigt model |4| (see Fig. bottom), with a power law 
spectrum of retardation times. The association seems 
natural, as it reflects the network, gelled, structure of 
the cytoskeleton, capable of supporting a certain amount 
of stress, without any full relaxation. It is known that we 
may mathematically translate the relaxation spectrum of 
the generalised Maxwell model to the retardation spec¬ 
trum of the Kelvin-Voigt model 130 ] [ 35 ]. The formu¬ 
las that allow us to make this passage are however non 
trivial, and the physical significance of the relaxation or 
retardation spectra may be lost. 

In this article, we will firstly highlight the fundamen¬ 
tal aspects of these interesting and universal weak power 
law behaviours, in the framework of the linear theory 
of viscoelasticity. This will help us to understand not 
only the reason for their universality, but also the reason 
for the large number of possible theoretical microscopic 
models that are able to explain them. Next, we will ex¬ 
tensively investigate two important types of microscopic 
models. The first, based on a generalised Maxwell model 
(see Fig. top), retains the fluid properties of the sys¬ 
tem and its relaxation modes. The second, based on a 
generalised Kelvin-Voigt model (see Fig. bottom), re¬ 
tains a more solid or gelled character, while presenting 
the weak power law behaviours. We will see that these 
two models, even though they have parallel descriptions, 
they translate essentially different microscopic physical 
systems. 

Using these models, we calculate with accuracy the 
power law exponents of G'{w) ^ and G"{w) ~ , 

associated with the power law spectra of relaxation and 
retardation times, % ~ l/r“ and C ~ r“, respectively. 
These depend on the exponent a and on the ratio be¬ 
tween the maximum and minimum relaxation or retar¬ 
dation times, Tmax and Tmin- Within the region cor¬ 
responding to 0 < a < 1, we find x > y ~ a and 
y > X ~ a for respectively the generalised Maxwell and 
Kelvin-Voigt models. The difference between the expo¬ 
nents X and y decreases as Tmin —> 0 and Tmax —>■ oo, but 
they are still substantial for the ratio Tmax/Tmin = 10® 


shown in Fig. and (to our knowledge, the accurate 
results shown in these figures were not described before). 
Beyond this region, for a < 0 or a > 2, we find very 
different and interesting exponents. 
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FIG. 1: Generalised Maxwell (top) and Kelvin-Voigt (bot¬ 
tom) models. 


POWER LAW RHEOLOGY 


In the linear theory of viscoelasticity [29], the relax¬ 
ation function G{t) is defined through the force response 
a{t) = eoG{t) to an applied step deformation £ = eoQ{t) 
(0(t) is the step function). Analogously, the creep com¬ 
pliance J{t) is defined through the deformation response 
e{t) = (JoJ{t) to an applied step force a = ao0(t). Also, 
the response of the system to an oscillatory deformation 
e = or force a = aoe™* is given, respectively, by 

a{t) = eoG*e^^* or e(t) = ctq where G* = G'+iG" 

is the complex modulus and J* = J' — iJ" is the complex 
compliance. G' (J') and G" (J") are defined as the elas¬ 
tic and viscous moduli (compliances). All these functions 
are related through: 


G*(w) = (iw)G(iw) = -=- 

{iw)J{iw) 


1 

J* (w) 


( 1 ) 


where the Laplace transform is given by /(s) = 

Jo From these relations it may be made clear 

that any relaxation function with a weak power law be¬ 
haviour: 

G{t) ~ ^ (0 < a < 1) (2) 

will necessarily lead to the weak power law behaviours 
of the elastic and viscous moduli. Indeed, the relaxation 
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function is a decreasing function of time, implying a > 0. 
For a < 1, its Laplace transform is given by G(s) ~ 
the complex modulus becomes G*(w) ~ {iw)°‘ = 
yielding 

G'{w) - G"{w) - ^ = tan ™ (3) 

It should be noted that, in these conditions (0 < a < 1), 
J(s) ~ l/s“+^ and J{t) ~ The system is considered 
to be solid like whenever 0 < a < 0.5 (as G" jG' < 1) 
and fluid like whenever 0.5 < a < 1 (as G"jG' > 1). 

In very general terms, any combination of iV viscoelas¬ 
tic elements could now be used as a basis for which a par¬ 
ticular distribution of viscoelastic parameters would lead 
to an approximation of the weak power law behaviour 
expressed in Eq. and consequently, of the weak power 
law behaviours of Eq. 

GENERALISED MAXWELL MODEL 

We consider first a generalised Maxwell (GM) model 
(see Fig. top) composed of a parallel array of N 
Maxwell viscoelastic elements. Each Maxwell element 
is composed of a dashpot of viscosity r\i in series with 
a spring of stiffness G^. It behaves elastically on small 
times scales and as a fluid on long ones. Its relaxation 
function is Gi(t) = Gie“*/'^b where = rji/Gi is the ele¬ 
ment relaxation time. Because the Maxwell elements are 
in parallel, the total relaxation function is just the sum: 

N 

GGM(i) = 5]G,e-*/"- (4) 

i=l 

This sum may be approximated by an integral by multi¬ 
plying it by dz = 1. If di is considered to be small, we 
have: 

/ ■’■JV JA 

(5) 

The derivative di/dr corresponds to the density of modes 
(per unit relaxation time). Choosing the change of vari¬ 
ables C, = I/r, and for t between the minimum and max¬ 
imum relaxation times, Tmin and r^ax, this integral ap¬ 
proaches a Laplace transform. To establish the weak 
power law behaviour G{t) ^ I/t“, the relaxation spec¬ 
trum must scale as 

= ( 6 ) 

Because we have two independent functions, Gr and 
di/dr, one of them may be chosen arbitrarily. 

An interesting and particularly simple choice consists 
of assuming a constant Gt = G. In this case, the condi¬ 
tion given by Eq. [^implies 

^ (C< = G) (7) 


Eor N sufficiently large, this spectrum will give the 
weak power law behaviours of Eq. G*{w) « 

G(zTir(;)“(ar(a)r(I—a)) (F is the gamma function), pro¬ 
vided 0 < a < I and 1 /ti = l/r^ax < w < I/T„i„ = 
l/(ri/iVi/“). 

This choice corresponds both to the Rouse and the 
Zimm models for unentangled polymers [3]. The Rouse 
model describes polymer chains through N Brownian 
beads in a viscous solvent (dashpot element) connected 
by springs. The chains have N vibrating modes with 
Ti = where Tfj is the Rouse relaxation time, 

V = 1/2 for ideal (random walk) chains and i/ « 3/5 
(Flory exponent) for self avoiding chains. The Zimm 
model extends the Rouse model by including hydrody¬ 
namic interactions mediated by the solvent. In this case, 
Ti = Tzi~^'^ where tz is the Zimm relaxation time. 
df = Ijv corresponds to the fractal dimension of the 
polymers. These models predict a = d//(2-|-d/) (a = 0.5 
for ideal chains) and a = df/3 respectively. 

The relaxation spectrum of Eq. also corresponds to 
the SGM model, which is based upon Bouchaud’s trap 
model of glassy dynamics [3T]. The essence of the SGM 
model may be described with simple arguments, as fol¬ 
lows. For each local element of the system, the energy 
landscape is composed of a large number of potential 
wells, or traps, with a density of deep traps that has an 
exponential tail, p{E) ^ (at the glass temper¬ 

ature, Xg = 1 in appropriate units). The equilibrium 
probability of finding at a certain moment an element 
in a trap of energy depth E is Peq{E) ^ t{E)p{E), 
where t{E) ~ is the average lifetime of the ele¬ 

ment, before it hops into a new trap, and a:: > I is a 
mean-field noise temperature (so Peg is no longer nor- 
malisable at the glass temperature). At characteristic 
times smaller than the average lifetime, the element be¬ 
haves essentially as a spring. For longer characteristic 
times, the element flows into the next trap. So, we may 
associate the behaviour of each element with a Maxwell 
element with a relaxation time that is equal to its av¬ 
erage lifetime. If we have a large number N of trapped 
elements, at a certain moment, the density of elements 
(per unit energy depth) is di/dE ^ Peq{E). The den¬ 
sity of elements (per unit relaxation time) is determined 
by di/dr = {di/dE){dE/dr) ~ t““. Taking the same 
quadratic shape for each energy well, G^ is considered to 
be constant, and the condition expressed in Eq. leads 
to a; = a+1. The system becomes solid like (0 < a < 0.5) 
as it approaches the glass transition. 

Had we chosen the power law behaviour Gi = Gi^ 
for the stiffenesses, the weak power law behaviours of 
Eq. could be recovered through the following discrete 
spectrum of relaxation times: 

’■* = = “b 

for 0 < a < I and l/r^nax < w < Train-, independently 
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FIG. 2: Elastic and viscous moduli, G'{w) (solid line) and 
G" (w) (dashed line) (in units of G) obtained from the gen¬ 
eralised Maxwell relaxation spectrum of Eq. (left) and 
the generalised Kelvin-Voigt retardation spectrum of Eq. 
(right), for a = —0.5,0.25,0.5,0.75,1.5 (top to bottom). 
Train = 10“®, Tmax = 10^ (iu arbitrary units), and N — 10^. 


of our choice of 13. As the continuous spectrum is even¬ 
tually the same as the one obtained for our simple previ¬ 
ous choice (corresponding to /3 = 0), the behaviour of all 
physical properties, in the large N limit, will be identical. 
Even the ranges beyond and r^ax will behave simi¬ 
larly, since in these limits, the elastic and viscous moduli 
depend essentially on one Maxwell element. Neverthe¬ 


less, the parameter (3 may play a role for the microscopic 
physical model. 

Mathematically, having one extra degree of freedom, 
the discrete relaxation spectrum of Eq. allows us to 
define not only the time relaxation limits, but also the 
number of elements we want to use, for any a. Indeed, for 
a given set of Tmax^ Train, ct, we are able to construct the 
weak power law behaviours of Eq. from a generalised 
Maxwell model with N modes by taking ri = Tmax and 

P = (Q;ln(Tmaa;/'rmm)/lniV) - 1. 

If we use the Laplace transform of Eq. we may 
determine without approximation the elastic and viscous 
moduli of the generalised Maxwell model, for any choice 
of discrete relaxation spectrum: 


N 

i=l 

N 

G'gm(^) = ^ Gi 

i=l 


{wTjf 
1 -I- {wTiY 

WTj 

1 -I- {wTiY ■ 


(9) 

( 10 ) 


By using these formulas, there are no restrictions on the 
parameter a. We have a well defined relaxation function 
(Eq. 1^ even for a < 0 or a > 1. 

Figure (left) shows the elastic and viscous moduli, 
G'qm{w) and for a = -0.5,0.25,0.5,0.75,1.5, 

Tmin = 10“^, Tmax = 10^ (in arbitrary units), and 
N = 10^. To obtain them, we have used the sums of Eq. 
I^andp^ together with the discrete relaxation spectrum 
defined by Eq. For 0 < a < 1, we recover the weak 
power law behaviours of Eq. Outside this range, we 
have for G' and G" different power laws. For 1 < a < 2, 
we have G' ^ w°‘ and G" ~ w^. For — 1 < a < 0, 
we have G' ^ and G" ^ w°‘ (this behaviour is ob¬ 
served near the glass transition of concentrated colloidal 
suspensions, and is associated with the slow relaxation 
modes caused by the cooperative dynamics of the open¬ 
ing of particle cages HO])- For Q;>2orQ!<—1, only one 
mode dominates, corresponding to the minimum and to 
the maximum relaxation times, respectively. These be¬ 
haviours had already been mentioned in the framework of 
the SGM model [53]. As these power law exponents (in¬ 
cluding those of Eq. correspond only to an approxima¬ 
tion (they would be valid for Tmin —> 0 and Tmax —t oo), 
we show in Fig. the numerically converged (we have 
increased N to verify it) computed values, established in 
the interior range 10“^ < w < 10^. It is important to 
mention that the exponent of G' is always greater than 
the exponent of G", even in the region 0 < a < 1 (see 
Fig. 1^, where they assume almost identical values. This 
is a distinctive feature of the generalised Maxwell model. 

It should be noted that we have performed the same 
calculations for different numbers of elements. The pos¬ 
sibility of changing /? allows us to have a good approxi¬ 
mation of the power law behaviours even for a reduced 
number of elements {N > 10). 
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FIG. 3: Exponents x and y vs a of the power laws G' ~ 
and G" ~ w^, obtained from the generalised Maxwell 
relaxation spectrum of Eq. |^(top) and the generalised Kelvin- 
Voigt retardation spectrum of Eq. [^(bottom), with Tmin = 
10“^, Tmax = 10^ (in arbitrary units), and N = 10®. These 
exponents were computed numerically, for every a, in the 
interior range 10~® < w < 10®. 


GENERALISED KELVIN-VOIGT MODEL 

We now turn to the generalised Kelvin-Voigt (GK) 
model (see Fig. bottom) composed of a serial ar¬ 
ray of N Kelvin-Voigt viscoelastic elements m- Each 
Kelvin-Voigt element is composed of a dashpot of vis¬ 
cosity rji in parallel with a spring of stiffness Gi. It be¬ 
haves elastically on long times scales and its dynamics 
comes from the viscous element. Its creep compliance is 
Ji{t) = Ji(l — e”*/®"’), where Ji = 1/Gi and = r]i/Gi is 
the element retardation time. Because the Kelvin-Voigt 
elements are in series, the total creep compliance is just 
the sum: 


N 

( 11 ) 

This sum may be approximated by the integral: 

Jgk- J J,(l-e-*/"0^dr (12) 

If we want to establish the weak power law dependence 
Jgk the retardation spectrum must then scale as 




(13) 


Again, because we have two independent functions, Jr = 
1/Gr and di/dr, one of them may be chosen arbitrarily. 
If we choose the power law behaviour Ji = Ji^ the 


condition given by Eq. 13 implies the following discrete 




FIG. 4: Exponents x and j/ vs a of the power laws G' ~ 
and G" ~ , obtained from the generalised Maxwell 

relaxation spectrum of Eq. [^(left) and the generalised Kelvin- 
Voigt retardation spectrum of Eq. (right). This is a blow 
up zoom of Eig. f or 0 < a < 1, where the weak power law 
behaviours of EqTl^ are approximately satisfied. 


spectrum of retardation times: 

Ti = (Ji = Ji^) 


(14) 


For N sufficiently large, this spectrum will give the 
weak power law behaviours of Eq. G*{w) « 

G{iTiw)°‘/{T{1 + a)F(l — a)), provided 0 < a < 1 and 
l/(riV(®+^)/“) = IjTmax < w < l/Tmin = I/ti- For a 
given set of Tmax, Tmin, C(, we are able to construct the 
weak power law behaviours of Eq. from a generalised 
Kelvin-Voigt model with N modes by making ti = Tmin 

and 13 = {a\n{Tmax/Tmin)/'^T^N) - 1. 

The exact elastic and viscous compliances are given 


through the Laplace transform of Eq. 11 


N 


^GKi'^) — ^ 


^ 1 -I- (wny 


N 


Jgk{'^) — T/ X 


WTi 


^ 1 -h {wnY 


(15) 

(16) 


Again, by using these formulas, there are no restrictions 
for the parameter a. In order to determine the elastic and 
viscous moduli we may simply use the relation between 
the complex modulus and compliance G*J* = 1. 

Figure (right) shows the elastic and viscous moduli, 
G'qj^(w) and G'l;j^{w) for a = -0.5,0.25,0.5,0.75,1.5, 
Tmin = 10“®, Tmax = 10® (in arbitrary units), and 
N = 10®. To obtain them, we have used the sums of 


Eq. ^ and together with the discrete retardation 
spectrum defined by Eq. [T^ For 0 < a < 1, we recover 
the weak power law behaviours of Eq. Outside this 
range, we have for G' and G” other power laws, which 
are essentially different from the power laws obtained for 
the generalised Maxwell model. For 1 < a < 2, we have 
G' ^ and G" ^ rc®. For — 1 < a < 0, we have 

G' ~ w® and G" ^ w~°' (this behaviour was observed in 
the simulations of a subcellular element model |32]). For 
a>2ora<—1, only one Kelvin-Voigt mode dominates. 
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corresponding to the maximum and to the minimum re¬ 
laxation times, respectively. As these power law expo¬ 
nents correspond only to an approximation, only valid as 
Tmin 0 and Tmax oo, we show in Fig. the nu¬ 
merically converged computed values, established in the 
interior range 10“^ < w < 10^. It is important to men¬ 
tion that the exponent of G" is always greater than the 
exponent of G", even in the region 0 < a < 1 (see Fig. 
1^, where they assume almost identical values. This is a 
distinctive feature of the generalised Kelvin-Voigt model. 
This inversion in the values of the slopes is reflected also 
in the intersection points of G' and G" (see Fig. right 
and left). If the intersection takes place at small values 
of w, for the generalised Maxwell model (when a < 0.5), 
it occurs at large values of w for the generalised Kelvin- 
Voigt model (for the same value of a) and vice versa. 

The generalised Kelvin-Voigt model, corresponding to 
an elastic system with well defined equilibrium lengths, 
seems more suitable to describe a gelled structure. To 
relate the parameter a with the physical properties of the 
microscopic system, we firstly assume that the density of 
elements with size I (per unit length) is given by the 
power law distribution di/dl ~ We also assume 

that the retardation time of an element of size I isr ^ P. 
So, the density of elements (per unit retardation time) is 
di/dr ~ Finally, if we take the general relation 

J , the condition of Eq. [Othat establishes the weak 
power law behaviours of Eq. [^implies a = {S — 

The parameter ^ has been related to the fractal dimen¬ 
sion of the gel, df [4]. This point is however subtle, and 
must be considered with care. Indeed, a gel with a well 
defined fractal dimension could be more easily identified 
with an irregular self-similar 3D network of Kelvin-Voigt 
elements, both in series and in parallel. The response 
of the Kelvin-Voigt elements along the direction of the 
oscillation will necessarily affect and be affected by the 
other 2 dimensions. In certain conditions, we were able 
to obtain the relation a = (<5 -|- df/3)/'y [ 55 ] . 

The viscosity of each element rji describes its resistance 
to the flow. It is reasonable to assume, as in Stokes’s law, 
that rji is proportional to the viscosity of the solvent and 
to the size of the element, r]i li. If we take Ji constant 
(as in the Rouse model, following from the equipartition 
theorem), ~ li, and we may obtain the interesting 
result a = df/3. There are other reasonable scaling laws. 
We may for instance invoke the idea of springs in series 
to justify Ji ^ li. In this case, Ti ^ If. 

CONCLUSION 

In this article, we have focused on the interesting weak 
power law behaviours G'{w) ^ G"{w) ~ (0 < a < 

1) that occur in a large variety of soft materials. By 
highlighting the relation between these power laws and 
a relaxation function that scales as G{t) ~ l/t“ (only 


when 0 < a < 1), we have stressed a fundamental reason 
for their universality. Any combination of N viscoelastic 
elements, with a particular distribution, or spectrum, of 
viscoelastic parameters may lead to these weak power law 
behaviours. 

They may appear in diluted, fluid systems, such as the 
generalised Maxwell model, that encompasses both the 
Brownian motion of elastic chains in a solvent, modelling 
unentangled polymers, or the important SGM model, 
that describes fluid behaviour above the glass transition. 
As the disorder and the metastability increase, the sys¬ 
tem becomes more solid like (whenever 0 < a < 0.5). 
The diluted bacterial cultures we have previously stud¬ 
ied seem to be satisfactorily described in this context. In¬ 
terestingly, in the range of angular frequencies in which 
we observed, in our cultures, the weak power law be¬ 
haviours, the exponent of the elastic modulus is always 
slightly greater than the exponent of the viscous modulus 

m- 

The weak power law behaviours may also occur in 
gelled, elastic systems, such as the generalised Kelvin- 
Voigt model. This model, although phenomenological, 
may retain two important characteristics of the cell’s cy- 
toskeleton, namely the idea of an existing pre-stress on 
its fibres, justifying the choice of using Kelvin-Voigt vis¬ 
coelastic elements which do not relax under stress, and 
the idea of a possible fractal structure, characterised by 
a power law distribution of elastic element lengths and 
retardation times [33] . Interestingly, some experimental 
observations seem to indicate a viscous modulus expo¬ 
nent greater than the elastic one [5]. The cytoskeleton 
is however a very complex active structure, and exhibits 
other important rheological properties that are not con¬ 
sidered here, such as stress stiffening [^ or a new power 
law exponent (of 3/4) for higher frequencies [35] . 

Particularly interesting, we find that beyond the re¬ 
gion corresponding to the weak power law behaviours 
(Q;<0ora>l), the generalised Maxwell and Kelvin- 
Voigt models give very different and distinctive rheo¬ 
logical behaviours, converging continuously to the 1- 
element Maxwell and Kelvin-Voigt models, respectively, 
for a < —1 or a > 2. 
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